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FIRST-ORDER EFFECTS OF PLANE SURFACES ON 


THE KINETIC BEHAVIOR OF GASES 

By Willard E. Meador 
Langley Research Center 

SUMMARY 

A linearized Boltzmann equation is solved exactly for the problem of constant veloc- 
ity and temperature gradients in a semi-infinite gas of Maxwellian particles. Contrary 
to similar previous work, the solution satisfies all physical requirements at the boundary 
without the introduction of complex exponential terms which are nonanalytic in the mean 
free path. The results apply also to simple Couette flow and are accurate over a range 
of Knudsen numbers extending well into the transition flow regime. Less rigorous exten- 
sions to large flow velocities are investigated in the midst of the gas and are compared 
with the work of Truesdell. 


INTRODUCTION 

The accommodation of gases to solid surfaces is a problem of long standing in 
kinetic theory. (See ref. 1 for a review of the earliest works.) Theoretical expressions 
for the associated velocity and temperature jumps in the immediate vicinity of a wall 
were first derived by Maxwell, but not without several assumptions which seemed to lack 
any rigorous support. In particular, the velocity distribution function was assumed to 
have the same characteristics close to the wall as in the midst of the gas, the velocity and 
temperature profiles were not permitted any unusual behavior, and the reflected particles 
were divided into two groups corresponding to specular and diffuse reflection. Only the 
last assumption is regarded as obsolete (see refs. 2 and 3) in the present paper. 

One of the questions concerning Maxwell's theory is whether the velocity and tem- 
perature jumps are real phenomena or simply mathematical consequences of the basic 
assumptions; for example, more rigorous treatments may yield continuous approaches of 
the flow velocity and gas temperature to the wall values and still satisfy the boundary con- 
ditions. Among the attempts to derive the behavior of gases near walls directly from 
kinetic theory is a series of three papers (refs. 4 to 6) by Wang Chang and Uhlenbeck 
which indeed show substantial trends toward the continuous approach to wall properties, 
but which also predict nonanalytic velocity distribution functions with respect to the mean 
free path. Such results are disturbing because they contradict both the Chapman-Enskog 



perturbation procedure (ref. 7) and the Grad moment expansion (ref. 8) in all but the near- 
continuum flow regime described by Navier-Stokes theory. Breakdowns of the standard 
techniques may not occur so rapidly with increasing values of the Knudsen number. 

A primary purpose of the present research is to show that Wang Chang and 
Uhlenbeck's nonanalytic distribution functions are direct consequences of the use of 
Maxwell's hypothesis about the character of the reflected particles. The combination of 
Boltzmann's kinetic equation with oversimplified boundary conditions inevitably leads to 
the incorrect appearance of nonanalytic hyperbolic terms. More physical boundary con- 
ditions with regard to the conservation of individual macroscopic moments, as opposed 
to Maxwell's separation of reflected particles into specular and diffuse groups, are shown 
in the present paper to be entirely consistent with the assumption that the velocity dis- 
tribution function has the same characteristics close to the wall as in the midst of the gas. 
Velocity and temperature jumps occur precisely as Maxwell predicted. 

These conclusions are derived from calculations based on the exact solution of a 
linearized Boltzmann equation for the problem of constant velocity and temperature gra- 
dients in a semi-infinite gas of Maxwellian particles. The linearization technique is 
equivalent to that of references 4 to 6; thus the only difference between the general 
methods is in the expression and use of constraints at the wall. Use also is made of the 
fundamental premise that exact solutions which satisfy all boundary conditions accurately 
describe the gas behavior. 

The present results apply also to simple Couette flow and are accurate for low 
Mach numbers over a range of Knudsen numbers extending well into the transition flow 
regime. This last feature provides a reliable extrapolation of the Navier-Stokes theory 
for the drag experienced by a flat plate; in addition, similar extrapolations are discussed 
for the heat transfer between flat plates at different temperatures. Finally, less rigorous 
extensions to greater flow velocities in the midst of a viscous gas are investigated and 
compared with Truesdell's work (ref. 9), which has been used in somewhat the same man- 
ner as that of Wang Chang and Uhlenbeck in criticizing the standard Chapman-Enskog and 
Grad approximations. (See ref. 10.) 


SYMBOLS 

a^^ distribution function constants 

A Wang Chang-Uhlenbeck accommodation coefficient 

B constant in Wang Chang-Uhlenbeck boundary condition 
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c particle velocity 

d parallel plate separation 

D drag force per unit area 

f velocity distribution function 

f D Maxwellian distribution function relative to mean gas velocity 

f o (0) Maxwellian distribution function with respect to gas temperature at surface 

of plate and zero flow 

G arbitrary velocity function 

h perturbation function relative to plate conditions 

AH excess random energy striking plate per unit area per second 

H 0 random energy striking plate per unit area per second 

H w random energy reemitted by plate per unit area per second 

A A 

i,j,k unit vectors in x-, y-, and z-directions 

k Boltzmann's constant 

K,KQ,Kj,Kg velocity parameters 
l mean free path 

m particle mass 

n particle number density 

nQ particle number density at surface of plate 


P 


scalar pressure 



p„„ third moment defined by equation (59) 

xzz 

JL 

p traceless pressure tensor 

q” conductive heat flux 

t time 

T temperature of gas 

Tq gas temperature at surface of plate 

T w plate (or wall) temperature if single plate 

Tj,T 2 temperatures of lower and upper plates 

if dimensionless particle velocity relative to mean gas flow 

U unit tensor 

v mean gas velocity 

v mean equilibrium particle speed 

v s slip velocity 

v w velocity of upper plate in Couette flow 

x,y,z Cartesian coordinates; designate vector and tensor components when used 

as subscripts 

a accommodation coefficient for energy 

/3 magnitude of dimensionless flow velocity 

/3 S dimensionless slip velocity 

dimensionless velocity (Mach number) of upper plate in Couette flow 
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y dimensionless particle velocity relative to stationary plate 

e temperature parameter defined by equation (30) 

e momentum parameter introduced in equation (24) 

N 

7] viscosity 

6 arbitrary velocity function 

A. thermal conductivity 

i; temperature parameter 

p mass density 

a accommodation coefficient for parallel momentum 

ct n accommodation coefficient for perpendicular momentum 

c t accommodation coefficient for parallel conductive heat flux 

c p perturbation function 

Special notations: 

collisional time derivative of velocity distribution function 
velocity average of G based on complete distribution function 
velocity average of G based on perturbation function f o (0)h 
velocity average of G based on Maxwellian distribution f o (0) 

Vector quantities without arrows denote magnitudes. 

THE LINEARIZED KINETIC THEORY 

The purpose of this section is to develop the linearized kinetic theory for a semi- 
infinite gas of Maxwellian particles having the velocity and temperature profiles 
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and 

T = T 0 (l + £z) (2) 

where K and £ are constants specified in the problem definition and Tg and /3 S 
are the gas temperature and dimensionless flow velocity at the surface (z = 0) of a sta- 
tionary flat plate. These expressions are assumed to hold throughout the gas. 

Equations (1) and (2) are employed directly in the Boltzmann kinetic equation 
(ref. 7) 

c • Vf = (1 + 0)c • Vf Q + f Q c • V0 = (3) 

for the steady-state velocity distribution function 

f = f 0 (l + 0) = ^ e" u2 (l + 0) (4) 

The symbol n refers to the particle number density and 



is the dimensionless particle velocity. Body forces such as gravity are neglected. 

If the macroscopic gas properties vary only with z, the gradient of the Maxwellian 
function f 0 can be written as 
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with the aid of equations (1) and (2) and the ideal gas lav/ p = nkT. The approximate 
equalities in equation (6) refer to the arbitrary neglect of quadratic and higher order 
terms in both /3 and (T - Tq)/t. 

This neglect of second and higher order terms, together with the assumption that 
the largest contributions to <p are first order (that is, linear in /3 or (T - To)/t), 
completely defines the linearization process and yields 


M r^/2r 2KUj 


Am ) 


^(u 2 _ 5\ + ldE + 8«£|u 

A 2/ p dz 9z z \9t/i 


( 7 ) 


from equations (3) and (6). The perpendicular distance from the plate is obviously lim- 
ited to moderate values by the restrictions on /3 and T. 

Another convenient assumption is that the diagonal elements of the traceless pres- 
sure tensor 

P"= 2p (<uu> - |<u 2 ) u) (8) 


are zero through first order, which is a familiar result in the midst of gases experiencing 
simple shear and implies a constant scalar pressure because of the steady-state equation 
of motion 



(9) 


The substitution into equation (7) of dp/dz = 0 and the trial function 


<p = ai u x u z + a 2 (u 2 - |)u z 


yields 


+ i(u 2 - |)lu z = (f) c = -| f o 


a l u x + - £ 


( 10 ) 


( 11 ) 


for Maxwellian particles (that is, interparticle potentials varying as the inverse fourth 
power of the separation). Only for such interactions do the Boltzmann collision integrals 
preserve the spherical harmonic character of equation (10) and permit an exact closed- 
form solution (ref. 7). The viscosity t? in equation (11) is a correct representation of 
the detailed collision dynamics. 

Equations (10) and (11) combine to give 


<P = - 




( 12 ) 


which is the exact first-order perturbation function in the midst of the gas and is assumed 
not to change in character as z approaches zero. However, this expression and the 
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velocity and temperature profiles of equations (1) and (2) also must satisfy the boundary 
conditions at the surface of the flat plate before they can be termed a precise physical 
description of the problem. Wang Chang and Uhlenbeck’ s constraints are not so satisfied. 
The critical subject of correct boundary conditions is discussed in some detail in the next 
two sections. 


BOUNDARY CONDITIONS 


Since the flat plate is stationary relative to the observer and located at the position 
where the gas temperature is Tg and the particle number density is ng, the most 
appropriate formulation of equations (4) and (12) for the specification of boundary condi- 
tions is 


x , /I . /n\/ m \3/2 / t oY / 2 f mc 2 / T l 

f - f 0 (l + 0) ~ n o(ng)^27TkTgj \T ) exp " 2kTg\T 

\3/2 / rnr»2'\ 


mc“ l •‘•CM , oof m \ 1 / 2 „ 

2kT 0 Uy + 2 ' 3 2kT 0 ) 1 




m 


1 + £&£? z + 2/3, 

2kTg P ^2kTg 


1/2 


(1 + 0 ) 


(1 + 0 ) 




1 + 2/3 


m_V 

2kT 0 j 


1/2 




'me' 


,2kTg 2 




(13) 


The expression for n employed in equation (13) is obtained from equation (2) and 
dp/dz = 0 to be 

n = n Q (l - £z) 

If the new definitions 

, m \V2_ 

y = f \ c 


2kTp 


and 


are introduced, and if the mean free path l is defined by the viscosity relation 
V = 


(14) 


(15) 


(16) 


(17) 


where p is the mass density and 


- 2 m) 1/2 

\ 7 ra 1 / 


( 18 ) 
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( 19 ) 


is the mean equilibrium particle speed, equation (13) becomes 
f = f o (0) (1 + h) 

with 


h = 2(0 S + Kz)y x + Jz(y2 - §) - + - §) 


jr z ( 20 ) 

The conservation laws to be satisfied by equations (19) and (20) at z = 0 are 



corresponding to particles, the x-, y-, and z -components of momentum, and energy, 
respectively, with gas-wall accommodation coefficients 0, a, a, a^, and a. Special 
definitions include AH for the magnitude of the excess random energy striking unit area 
of the wall in 1 second, the parameter relating to the z-momentum associated with 

AH, and the symbolic integral 

<G|h> = A- J f o (0) Gh dc = -A_ J e -^ 2 Gh df (26) 

for velocity-averaged quantities. 

If each accommodation coefficient refers on the average to the fraction of a parti- 
cle's corresponding macroscopic property which is lost on impact with the wall, the left- 
hand sides of equations (21) to (25) represent the total amounts of such properties absorbed 
in 1 second by a unit area of the wall. The steady-state conservation laws require these 
absorptions to be balanced by the net property fluxes on the right-hand sides of equa- 
tions (21) to (25). Only in equation (24) and AH in equation (25) remain to be speci- 

fied in order to utilize these conditions. 
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Since the random energy of the incident stream is governed entirely by the 
Maxwellian function f o (0) at z = 0, the magnitude of such energy crossing unit area 
in 1 second must satisfy 


tt _ m 

n °~T m 


s 

o\- 


2kT n J y>° 


-y 2 o 

e y*y z dy 


L 


r z >o 


-y2 

; y z dy 


n 0 /2kT 0 \ 1/2 


71-3/2 \ m 


f e ~ y 2 y z d y 
J r z > o z _ 


= ( 2kT o)hf) (27) 


The terms in the two parentheses following the second equality sign in equation (27) cor- 
respond, respectively, to the average energy per particle issuing toward the wall from an 
equilibrium gas at temperature Tq and the number of such particles per unit area per 
second. A similar reflected stream at the temperature T w of the wall (but with the 
same number of equilibrium particles per unit area per second in order not to imply a 
net particle flux (see ref. 1, p. 313)) would carry energy away in the amount 

H„ = (28) 

The substitution into equation (25) of the energy difference 


AH = H 0 



(29) 


obtained from equations (2), (27), and (28), the definition 

e = T 0 " T w (30) 

T 0 

and the identity 


^ r rz<0 e " y2y2rz dy " ~"~ x/2 <si) 

yields the final energy condition 



This particular form is chosen for subsequent comparisons with Wang Chang and 
Uhlenbeck’s constraints. 

Of the five boundary conditions listed in equations (21) to (24) and (32), only the con- 
servation of energy and the x-component of momentum are not satisfied automatically by 
the h of equation (20) if 
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Equation (33) is obtained by using equation (20) to evaluate the integrals in equation (24) 
and solving the resulting expression for e^. The proportionality to | confirms the 
interpretation that e-^ is related to AH and vanishes in the absence of temperature 
gradients. 

Accordingly, equations (22) and (32) are sufficient to determine the two remaining 
unknown parameters j3 s and e and thereby complete the first-order description of a 
semi-infinite gas with given velocity and temperature gradients and mean free path and 
with given wall and gas -wall parameters T w > o', and o^. Each given quantity either 
contributes to the basic problem definition or must be determined from the details of gas- 
wall interactions; hence, the role of pure kinetic theory is finished with the application of 
equations (22) and (32) to equation (20). 

Pertinent integrations for this purpose are given by the expressions 

(34) 


( y x y z| h ) ~ 71-1/2 
( y2y z| h ) yz<0 = 

( y2y z| 1 ) yz < 0 = “ 7r “ 1/2 

and 

the appropriate combinations of which are 
§(0s + Kz) = K2 

and 

according to equations (22) and (32). 

Hence, the velocity and temperature jumps, v g and AT, satisfy 

v _ /^kT^/Z _ (2 - a)l(2 kT^ 1 / 2 ^ _ (2 - o)Z dv 

s \m/ a \ m / K a dz 


(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 
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and 


AT = T 0 - T w = T 0 e 


15(2 - a)g T q _ 15(2 - a )i dT 
8a 8a dz 


with the aid of equations (1), (2), and (30). 


(42) 


These results are precisely those of Maxwell (see ref. 1, pp. 296 and 314) but with 
the distinction that they are now supported by a rigorous first-order kinetic theory of 
Maxwellian particles. This conclusion depends, of course, on the completeness of the 
boundary conditions of equations (21) to (24) and (32). The question might be asked, for 
example, whether other particle properties should be constructed from the basic ones of 
particle identity, momentum, and energy, and then forced by steady-state arguments to 
satisfy their own conservation laws at the wall. 


A typical such constructed property is the x-component of a particle's heat flux 
which in the midst of a gas seems to infer a net flow of this property in the z-direction 
according to the integral ^^x^z |^)‘ a new accommodation coefficient a' be 

introduced so that 


a' 



<r 2 y x r z 


l h > 


(43) 


represents an additional constraint on the solution of the kinetic equation at z = 0? Or 
should a' be defined exclusively by equation (43) without any independent existence from 
detailed gas-wall physics? The most likely interpretation is the following: Since the 
conservation of energy and momentum completely describes the dynamics of gas-wall col- 
lisions if the accommodation coefficients a, a, and are known, equation (43) serves 

only to express o' in terms of a, a, and o-kt and thus should not be regarded as an 
additional restriction on the velocity distribution function. More precisely, 


a' = 


14a 
a + 12 


(44) 


from the application of equation (20) to equation (43) and the subsequent use of equa- 
tions (1) and (41). 

Although this distinction between the conservation of basic and constructed particle 
properties may seem trivial from the standpoint of a straightforward treatment of the 
present problem, it is absolutely essential for the critical analysis of Wang Chang and 
Uhlenbeck’s boundary conditions given in the next section. Their application of equa- 
tion (43) with a' equal to a instead of the value in equation (44) is directly responsible 
for the previously mentioned nonanalytic exponential contributions to h in references 5 
and 6. 
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I 


WANG CHANG AND UHLENBECK'S THEORY 

Wang Chang and Uhlenbeck* s theory (refs. 4 to 6) of simple shear in a semi-infinite 
gas of Maxwellian particles differs from that of the preceding sections only in the replace- 
ment of the conservation laws of equations (21) to (24) and (32) by the single microscopic 
condition 

h(y z ) = AB + (1 - A)h(-y z ) (y z > 0) (45) 

at z = 0. The accommodation coefficient A is said to represent the fraction of parti- 
cles striking the plate which is reemitted with the Maxwellian distribution corresponding 
to Tw and zero flow velocity - that is, which is deprived by the plate of all properties 
characteristic of the incident stream. The remaining particles are specularly reflected, 
whereas the unknown constant B is related to the gas density near the plate. These 
concepts are reiterated in a later survey by Uhlenbeck and Ford (ref. 11). 

But equation (45) must be statistically equivalent to the definition of A as the aver- 
age fraction of a particle's total nonequilibrium properties (plus excess random energy 
and associated z -momentum) absorbed by the plate; hence, the question must be asked 
whether a single accommodation coefficient is meaningful. In particular, the 
e(y 2 ,y x ,y y) y z - moment equation (45) yields 

^(y 2 ,y :i ,7'yjy 2 |li),, z>0 = AB <^(r 2 >V y x M}> rz >o + (1 ' A) 4 y ' 2 ^ x ^ yz ! h '' Tz Wz >0 

= -AB<e( y 2,r x ,r y )r z |l) yz<0 + (A - i><e(r 2 ,r x ,r y )r z |h) Vz<0 (46) 


which combines with the identity 



The coefficient A clearly refers in this expression to the average fraction of a 
particle's #(y2,y x ,yy)-property absorbed by the plate and h is the perturbation function 
of Wang Chang and Uhlenbeck. But equation (48) also raises the obvious and critical 
question of whether Wang Chang and Uhlenbeck' s constants A and B are adequate to 
describe the infinite variety of absorptions associated with different choices of 6. 
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Similarly, the y^-moment of equation (45) gives 

A ^Z 2 | h /y z <0 " B ( y Z 2 l 1 )y z <o) = {^Z 2 | h )y z <0 " {^Z 2 l^y z >0 ( 49 ) 

for the conservation of the z -component of momentum. Again the dependence of A on 
the moment taken is obvious: Unless the same average absorption fraction applies to 
every particle property, Wang Chang and Uhlenbeck's microscopic boundary condition is 
incorrect. 

However, the error is not as serious in some applications as might be inferred 
from the preceding discussion. A detailed comparison between equations (48) and (49) 
and equations (21) to (24) and (32) is given in table I for the present problem and also for 
the problems of zero shear (dv/dz = 0) and no temperature gradient. The h of equa- 
tion (20) is used in the calculations of A and B, which thus correspond to the exact 
first-order solution of the present research. 


TABLE I.- CORRELATIONS OF WANG CHANG AND UHLENBECK'S 
BOUNDARY CONDITIONS WITH THE CONSERVATION LAWS 
FOR THREE DIFFERENT PROBLEMS 


Property 

conserved 

_ 

dz 

— = 0 
dz 

dz 

0 

dz 

dz ’ 

dz 

A 

B 

A 

B 

A 

B 

Particles 

Arbitrary 

0 

Arbitrary 

0 

Arbitrary 

0 

x-momentum 

a 

Arbitrary 

Arbitrary 

Arbitrary 

a 

Arbitrary 

y-momentum 

Arbitrary 

Arbitrary 

Arbitrary 

Arbitrary 

Arbitrary 

Arbitrary 

z -momentum 

Arbitrary 

0 

°N 

" e N 

°N 

i “ e N 

Energy 

Arbitrary 

0 

a 

-e 

a 

-6 


The most outstanding feature of table I is the liberal sprinkling of arbitrary values 
of A and B for each problem. Although two of the problems require more than one 
value of A and B and thereby invalidate Wang Chang and Uhlenbeck's boundary condi- 
tion of equation (45), the one corresponding to simple shear with no temperature gradient 
clearly shows the constants A = a and B = 0 to be consistent with all five conservation 
laws. Why then do not Wang Chang and Uhlenbeck (refs. 5 and 6) obtain the results of this 
paper when T is constant? What is the origin of their nonanalytic exponential contribu- 
tions to the velocity distribution function and the velocity profile? 


14 



Answers to these questions are found in Wang Chang and Uhlenbeck’s special 
applications of equation (45) rather than in the properties listed in table I. Just as in the 
preceding sections, the boundary conditions of references 4 to 6 are applied only in the 
integrated forms of equations (48) and (49), but with no restrictions on the function 
q(q2 ,r x ,y y ). Consequently, to take one example, the function d = y^y x is both permitted 
and actually used by Wang Chang and Uhlenbeck to obtain the constraint 

cr ( y2y x^z| 1 )y z <0 = <r 2 y x r z |h) (50) 

for A = ct and B = 0 in equation (48). Except for the nonphysical values cr = 0 and 
ct= 2, this expression is clearly inconsistent with equations (43) and (44) and thus imposes 
an improper condition on the solution of the kinetic equation. The single accommodation 
coefficient A in equation (45) is never sufficient to describe the absorptions by the wall 
of all gas properties or conceivable moments. 

One of the principal objections to standard kinetic theories is removed by this anal- 
ysis of Wang Chang and Uhlenbeck' s research. No longer can the basic assumptions 
underlying the Burnett or 13-moment formulations be dismissed on the grounds that they 
disagree with the first-order behavior of a viscous gas near a wall. (See ref. 10.) Other 
difficulties do exist (ref. 12) and some criticisms are justified, especially with regard to 
the higher Chapman-Enskog approximations, but such aspects are irrelevant to the simple 
problems considered here. The exact first-order distribution function of equations (4) 
and (12) is precisely Grad's 13-moment solution. 

More recent work in the same problem area includes two papers by Scharf (refs. 13 
and 14), in which boundaries are regarded as sources of energy and momentum in the 
gaseous kinetic equations. The source terms appear as combinations of Dirac delta func- 
tions centered on the walls and correspond to gas properties at such positions having 
average values between those of the wall and those of the gas in the absence of the wall. 
These additions of explicit mixtures of gas and wall properties to the gaseous kinetic 
equations are neither required nor implied by kinetic theory and boundary conditions. 

COUETTE FLOW 


The results of the preceding sections are easily extended to the problem of simple 
Couette flow between two parallel plates, one of which is stationary at z = 0 while the 
other at z = d moves in the positive x-direction with the constant velocity v w . Equa- 
tion (41) and the linear velocity profile 


v = v s + z 


dv 

dz 


z + - 


(2 - o)l 


a 


(51) 
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are rigorously valid through first order, so that 


dv _ ( v w ~ v s) ~ v s _ v w 
dz d d 

Finally, since the traceless pressure tensor element P xz is -t) dv/dz from 
equations (8) and (12), the drag force D exerted on unit area of the lower plate must 
satisfy 


1 + 


2(2 - o H 
od 


( 52 ) 


D = -°P XZ = ^ 


i , 2(2 - a)f 


od 


according to equation (52). 


(53) 


If the temperature gradient is zero, the only restriction imposed on equations (51) 
to (53) is that the Mach number /3 W of the upper plate be small - or more precisely, 
from equation (12) , that the parameter Z/3 w /d be small. Hence, the Knudsen number 
l/d can extend almost to unity in many applications, which is well above the lower limit 
usually assigned to the transition flow regime, and equation (53) represents a correct 
generalization of the Navier-Stokes formula 


D = 


^w 

d 


(54) 


A similar study of the conductive heat flux between two stationary parallel plates, 
one at z = 0 with temperature Tj and the other at z = d with temperature T£, shows 


= gMV/2/2-x „ dT . dT = .j X ( T 2- T l) 

H ^\m ) \ / 4m dz dz d 


1 + 


15(2 - g)i 
4ad 


1-1 


(55) 


to be restricted only by small values of Z(T2 - Tjj/Td. The successive steps in equa- 
tion (55) are obtained by using equations (2), (12), and (42). 


HIGHER MACH NUMBERS 


In addition to the criticism based on Wang Chang and Uhlenbeck's theory, Schaaf 
and Chambre (ref. 10) state that the Burnett and 13-moment representations do not ade- 
quately approximate the results found by Truesdell (ref. 9) for the less difficult problem 
of simple shear in an infinite medium of Maxwellian particles. Since Truesdell claims 
an exact solution of the complete Boltzmann equation, at least a partial analysis of his 
work is appropriate to the present research. 

Truesdell’ s method begins with the complete set of macroscopic equations of change 
obtained by taking all moments of the rigorous Boltzmann kinetic equation. However, 
instead of following Grad's example of using an approximate distribution function to 
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break the chain of relations after a specified number of moments, he accomplishes the 
same purpose by making arbitrary assumptions about the spatial or time variations of all 
but a few of the infinite number of moments. In the particular problem of simple shear, 
Truesdell assumes all second and higher moments to be spatially homogeneous; as a 
result, only the equations of continuity, motion, and energy balance are nontrivial. The 
additional assumptions of constant number density and a linear velocity profile then force 
the problem to be time dependent, a fact which Truesdell uses to criticize previous 
steady-state solutions based on Burnett's equations. 

The contention of the present paper is that Truesdell' s criticism is invalid for two 
reasons: first, the higher moments are not subject to direct laboratory control (or 
accurate theoretical speculation) and may behave very differently from his assumptions; 
secondly, Truesdell's criticism fails to distinguish between the problem conditions arbi- 
trarily imposed by his assumptions and the somewhat different Conditions underlying the 
Burnett solutions. If, for example, the simple shear is maintained by the relative motion 
of two widely separated parallel plates to establish a definite plane of symmetry, the gas 
temperature will eventually peak at this plane and cause a steady-state conduction of 
friction-generated heat toward each plate. Truesdell's restriction to spatially homoge- 
neous second moments (including the temperature) precludes this heat conduction and thus 
corresponds to a steadily increasing temperature with time. 

A better and certainly more physical approach is to insist upon a steady state and 
deduce from the kinetic theory what the velocity profile must be. If the expansion param- 
eter l d/3/dz is small enough, the linearized theory of the preceding sections is appli- 
cable and predicts a constant velocity gradient. For somewhat larger values of l d/3/dz 
(or the maximum Mach number if the flow is Couette), a convenient procedure starts with 
the following macroscopic equations for spatial variations in the z-direction and a mean 
particle velocity in the x-direction: 

Conservation of x-momentum: 

£> X z = Constant (56) 


Conservation of z-momentum: 


P zz + p = Constant 
o 

and the P xz equation: 


o 

p xz + 




(57) 


(58) 
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where 


^xzz 



(59) 


Each of equations (56) to (58) is obtained from a moment of the Boltzmann equation. 


Unless higher moment equations are considered, some assumption must be made 
about the character of p xzz before equation (58) can be utilized effectively. If the sub- 
sequent analysis is restricted to third order, which at most is the order of p in this 
problem, and if the leading contribution to any gas property is assumed to be constant or 
proportional to z, then P xzz satisfies 


d^p 

*xzz 

dz^ 


= 0 


(60) 


Although the assumption underlying this restriction is far from established, it at least 
conforms to the previous first-order results for the flow velocity, pressure tensor, and 
heat flux. 

The second and final assumption is that the third-order flow velocity must satisfy 


v = Kq + Kjz + I^zS 


(61) 


where Kj and K 3 are first- and third-order constants, respectively. 

Since Kj is input data controlled by whatever mechanism maintains the flow, 
the determination of K 3 as a function of Kj will complete the present description. 
Accordingly, the differentiation of equation (58) and use of equations (56), (57), (60), 
and (61) yield the third-order expression 

U ^)- 6 ” KiZ+K ^ =0 m 


and thus 


k 3 = - 
6 6rjz dz 


(63) 


o 

for constant p and second- order P zz . The coefficient K 3 is obviously a third-order 
constant because 77 is proportional to some power of the temperature and the tempera- 
ture undergoes a second-order parabolic decay from its peak value at the symmetry plane. 
No temperature gradient was induced at the first-order level in the previous sections. 


The two fundamental conclusions of this brief study are summarized as follows: 

( 1 ) induced temperature gradients cause spatially dependent viscosities which lead to 
modifications of the velocity profiles in the simplest kinds of viscous flow, and ( 2 ) these 
modifications are at most third order in the expansion parameter £ or l d/3/dz. Hence, 
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the Burnett prediction of simple shear (that is, a linear velocity profile) in a steady- 
state infinite medium is entirely consistent with the second-order limitations of that 
approximation. 


CONCLUDING REMARKS 

The simple problems of viscous flow near a flat plate, Couette flow between parallel 
plates, and heat transfer between parallel plates are solved exactly through linear terms 
in typical Chapman-Enskog expansion parameters. Accurate corrections for Knudsen 
numbers extending almost to unity are applied to the Navier-Stokes theory of Couette flow 
with low Mach numbers. None of these examples displays the nonanalytic exponential 
character of velocity distribution functions and calculated results found near the walls by 
Wang Chang and Uhlenbeck. Differences between the two theories are explained in terms 
of boundary conditions, some of which are nonphysical in Wang Chang and Uhlenbeck' s 
procedure. Less rigorous extensions of the present theory to higher Mach numbers show 
that Truesdell's disagreements with the standard Chapman-Enskog and Grad approxima- 
tions are caused by nonphysical assumptions about the behavior of higher moments. 

In summary, much of Maxwell's original theory of velocity and temperature jumps 
is supported by the present exact solutions of linearized Boltzmann equations for 
Maxwellian molecules; more specifically, his assumptions about the distribution function 
and the nature of velocity and temperature profiles in the immediate vicinity of a wall are 
accurate through first order. Only the hypothesis that the reflected particles can be 
divided into diffuse and specular groups is obsolete. 

Extensions to more general interparticle interaction potentials and system geome- 
tries can be very difficult, but no major changes are expected in the basic conclusions. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 25, 1970. 
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